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We compute the spectrum of primordial gravitational wave perturbations 
in open de Sitter spacetime. The background spacetime is taken to be the 
continuation of an 0(5) symmetric instanton saddle point of the Euclidean 
^ no boundary path integral. The two-point tensor fluctuations are computed 



directly from the Euclidean path integral. The Euclidean correlator is then 
analytically continued into the Lorentzian region where it describes the quan- 



turn mechanical vacuum fluctuations of the graviton field. Unlike the results 



of earlier work, the correlator is shown to be unique and well behaved in the 
infrared. We show that the infrared divergence found in previous calculations 
is due to the contribution of a discrete gauge mode inadvertently included in 



Q_i' the spectrum. 

43 



I. INTRODUCTION 

One appeal of inflationary cosmology is its mechanism for the origin of cosmological 
perturbations. The de Sitter phase of exponentially-rapid expansion quickly redshifts away 
any local perturbations, leaving behind only the quantummechanical vacuum fluctuations in 
the various fields. During inflation, these perturbations are stretched to macroscopic length 
scales and subsequently amplified, to later seed the growth of the large scale structures in 
the present-day universe. A particularly clean example of this effect are the gravitational 
wave perturbations of the spacetime itself. These tensor perturbations contribute to the 
cosmic microwave background anisotropy via the Sachs- Wolfe effect. They may potentially 
provide an observational discriminant between different theories of open (or closed) inflation 
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because their long-wavelength modes strongly depend on the boundary conditions at the 
instanton that describes the beginning of the inflationary universe [|l|. 

Although the tensor spectrum has been successfully computed in realistic 0(3, 1) invari- 
ant models for an open inflationary universe fl[], the problem of calculating the primordial 
gravitational waves in perfect open de Sitter spacetime has remained a paradox for some 
time. The previous literature claims that the spectrum of gravitational waves in perfect de 
Sitter space is infrared divergent for all physically well-motivated initial quantum states of 
an eternally inflating universe @-f|. Breaking the 0(4, 1) invariance of de Sitter space by 
going to a realistic inflationary model introduces a potential barrier for the tensor fluctu- 
ation modes, and it has been argued that the bubble wall acts to regularise the divergent 
spectrum in perfect de Sitter space ||. 

Previous calculations of the gravitational wave spectrum in open de Sitter space 
are based on a mode-by-mode analysis. One has a prescription for the vacuum state of the 
graviton that is imposed on every mode separately, on some Cauchy surface for the de Sitter 
spacetime. Then one propagates each mode into the open universe region. In this paper we 
instead compute the two-point tensor correlator in real space. In doing so, we have obtained 
an infrared finite tensor spectrum. The difference in the two approaches is related to the 
non-uniqueness of the mode decomposition in an open universe, as we shall explain. 

As an aside, we mention in this context that also fluctuations of a massless minimally 
coupled scalar field in de Sitter space do not break 0(4,1). In some prior literature (see 
e.g. fLlf ) it is shown that there is no de Sitter invariant propagator for such a scalar field. 



However, the scalar field is not itself an observable since the action depends only on its 
derivative, and there is a symmetry — > <p+ constant. In fact, correlators of space or time 
derivatives of <j) are de Sitter invariant, and since these are the only physical correlators in 
the theory, de Sitter invariance is unbroken. 

We implement the Hartle-Hawking no boundary proposal || in our work by 'rounding 
off' open de Sitter space on a compact Euclidean instanton, namely a round four sphere. The 
fluctuations are computed in the Euclidean region directly from the Euclidean path integral, 
to first order in h around the instanton saddle point. The Euclidean two-point correlator is 
analytically continued into the Lorentzian region where it describes the quantum mechanical 
vacuum fluctuations of the graviton field in the state described by the no boundary proposal 
initial conditions. There is no ambiguity in the choice of initial conditions because the 
Euclidean correlator is unique. 



II. TENSOR FLUCTUATIONS ABOUT COSMOLOGICAL INSTANTONS 

In quantum cosmology the basic object is the wavefunctional *ff[hij, <j)], the amplitude for 
a three-geometry with metric hij and field configuration <fi. It is formally given by a path 
integral 

f 13, \vg}[V<P}e^\ (1) 

Following Hartle and Hawking || the lower limit of the path integral is defined by 
continuing to Euclidean time and integrating over all compact Riemannian metrics g and 
field configurations <fi. If one can find a saddle point of ([!]), namely a classical solution 
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satisfying the Euclidean no boundary condition, one can in principle at least compute the 
path integral as a perturbative expansion to any desired power in %. 

In this paper we wish to compute the two-point tensor fluctuation correlator in open de 
Sitter spacetime, 

ds 2 = -dt 2 + sinh 2 (t) (d X 2 + smh 2 { X )dQ 2 2 )) . (2) 

Open de Sitter space may be obtained by analytic continuation of an 0(5) invariant instan- 
ton, describing the beginning of a semi-eternally inflating universe. The analytic continu- 
ation is given by setting t = —io and the radial coordinate \ = iQ, where Q is the polar 
angle on the three sphere (see M). The instanton obtained in this way is a solution of the 
Euclidean equations of motion with the maximal symmetry allowed in four dimensions. It 
takes the form of a round four sphere with line element ds 2 = da 2 + sm 2 (a)dQ 2 , where dVL\ 
is the line element on S 3 . It is useful to introduce a conformal spatial coordinate X defined 
by f^ 2 so that the line element takes the form 

J J cr sin ct' ' 

ds 2 = cosh" 2 X (dX 2 + dnfj . (3) 

On the four sphere X then ranges from — oo to +oo. 

The principles of our method to calculate cosmological perturbations are described in 
detail in Jl]||. The instanton solution provides the classical background with respect to 
which the quantum fluctuations are defined. In the Euclidean region the exponent iS in 
the path integral becomes —Se = — (So + S2), where Se is the Euclidean action, So is 
the instanton action and S2 the action for fluctuations. We keep the latter only to second 
order. The path integral for the two-point tensor fluctuation about a particular instanton 
background is then given by 

ff Mf ry u _ WM^^Mj (a) 

[tijWi'A* )) - J[V5g][V5<l>]e-* ' [) 

To first order in h the quantum fluctuations are specified by a Gaussian integral. The 
Euclidean action determines the allowed perturbation modes because divergent modes are 
suppressed in the path integral. The Euclidean two-point tensor correlator is then ana- 
lytically continued into the Lorentzian region where it describes the quantum mechanical 
vacuum fluctuations of the graviton field in the state described by the no boundary proposal 
initial conditions. 

To find the perturbed action S2 that enters in the path integral (f|), we write the per- 
turbed line element in open de Sitter space as 



ds 2 = sinh- 2 (r) (-(1 + 2A)dr 2 + Stdx'dr + ( 7ij + h^dx'dx 3 ) , (5) 

where the fields A, Si and hy are small perturbations. Because we are interested in the 
gravitational wave spectrum in the open slicing of de Sitter space, we will only retain 0(3, 1) 
invariance in our calculation. 

The quantities Si and may be uniquely decomposed as follows |TTJ, 



hij = \hli 3 + 2 (ViVj - ^A 3 ) E + 2F m + Uj, 
Si = B\ t + Vi. (6) 
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Here A3 is the Laplacian on S 3 and \j the covariant derivative on the three-sphere. With 
respect to reparametrisations of the three-sphere, h, B and E are scalars, Vi and Fi are 
divergenceless vectors and is a transverse traceless symmetric tensor, describing the grav- 
itational waves. Because gauge transformations are scalar or vector, the perturbations 
are automatically gauge invariant. 

It is important to note that the gauge invariance of follows from the uniqueness of 
the above decomposition. This is only true however for bounded (asymptotically decaying) 
perturbations |10| . If one does not impose suitable asymptotic conditions on the fields, a 



degeneracy appears between scalar and tensor perturbations that introduces a discrete gauge 
mode in the tensor spectrum, which plays a crucial role in the divergent behaviour of the 
correlator. We come back to this point in Section V. 

We now substitute the decomposition (|]) into the Lorentzian action for gravity plus a 
cosmological constant, 

S=±-f d A x^9 (R - 2A) - - / d 3 x^K, (7) 

The scalar, vector and tensor quantities decouple. Keeping all terms to second order, we 
continue the perturbed Lorentzian action to the Euclidean region. The scalar and vector 
fluctuations are pure gauge in perfect de Sitter space. The tensor perturbations £y yield the 



following well-known positive Euclidean action ||12|| : 

1 r^vi 



Here prime denotes differentiation with respect to the conformal coordinate X. After per- 
forming the rescaling = — %-y and integrating by parts we obtain 



S 2 = — f d 4 x^t l3 (K + 3- A 3 ) t ij + — 17 d 3 x^t i:j t ij tanh(X) 

OK, J ^ ' oK IJ 



(9) 



where the Schrodinger operator 



k -^-^rr-^ +ui - x) - < 10) 

Because the fluctuations are specified by a Gaussian integral, we can solve the path in- 
tegral (Q) by looking for the Green function of the operator in its exponent. The potential 
U(X) for the fluctuation modes is well known to be perfectly reflectionless. However, chang- 
ing its shape slightly would introduce some reflection which becomes increasingly significant 
at small momenta. Such a change corresponds to breaking the 0(5) invariance of Euclidean 
de Sitter space and is exactly what happens in the 0(4) invariant Hawking-Turok and 
Coleman-De Luccia instantons that describe the beginning of realistic open inflationary 
universes. This difference between both classes of instantons has profound implications for 
the tensor perturbations about them, especially for their long- wavelength regime 0. The 
operator K has in all three cases a positive continuum starting at eigenvalue p 2 = 0, as well 
as a single bound state = b(X)qij(Q) at p = i which turns out to be a trivial gauge mode. 
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III. THE EUCLIDEAN GREEN FUNCTION 



To evaluate the path integral we first look for the Green function G^ % 3 (X, X', f2, 
of the operator in (0). The Euclidean fluctuation correlator (|j) will then be given by 
cosh(X) cosh(X')G^ 31 3 . The Euclidean Green function satisfies 

^ [k + 3 - a 3 ) Gg Vj ,{x, x', n, n') = s(x - x^-W^n - n'). (11) 

If we think of the scalar product as defined by integration over S 3 and summation over 
tensor indices, then the right hand side is the normalised projection operator onto transverse 
traceless tensors on S 3 . 

The Green function G% \,y can only be a function of the geodesic distance fl') if it 
is to be invariant under isometries of the three-sphere. This suggests that 

+JOO 

G«<, ,.,(/*, X, X>) = 4/c £ G p (X, X')W { ;i ,.,(/.), (12) 

p=3i 

where (//) is a bitensor that is invariant under the isometry group 0(4). It equals 

the sum ( |A2p of the normalised rank-two tensor eigenmodes with eigenvalue X p = p 2 + 3 of 
the Laplacian on S 3 . Note that the indices i,j lie in the tangent space over the point Q 
while the indices i',f lie in the tangent space over the point Q'. On S 3 we have 

The motivation for the unusual labelling of the eigenvalues of the Laplacian is that, as 
demonstrated in the Appendix, in terms of the label p the bitensor on S 3 has precisely the 
same formal expression as the corresponding bitensor on H 3 . It is precisely this property 
that will enable us in Section IV to continue the Green function from the Euclidean instanton 
into open de Sitter space without decomposing it in Fourier modes. The relation between 
the bitensors on S* 3 and H 3 together with some useful formulae and properties of maximally 
symmetric bitensors are given in Appendix A. 

Since the tensor eigenmodes of the Laplacian on S 3 form a complete basis, we can also 
write 

+ioo 

rh%, 3> {n - no = E KU>M^'))- ( 14 ) 

p=3i 

Hence by substituting our ansatz ( |i~2|) for the Green function into (0) we obtain an equation 
for the X-dependent part of the Green function, 

(K - p 2 ) G P (X, X') = 5(X - X'). (15) 
The solution to equation fll^D is 



g p (x, x') = ±- \%(x)¥ p (x')G(x - x') + ¥ p (x)%(x')e(x' - x)} . (ie) 



A p 
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^fp(X) is the solution to the Schrodinger equation that tends to e~~ tpX as X — > — oo, and 
typ(X) is the solution going as e ipX as X — > +00. The factor A p is the Wronskian of the 
two solutions. Since the potential is reflectionless on the round four sphere the left- and 
right-moving waves do not mix and they equal the Jost functions g± p {X) with nice analytic 
properties. The solutions may be found explicitely and are given by 

f %(X) = (tanhX-^)e ipX 

{ W p (X) = (t&nhX + ip)e~ ipX 1 ' 

and their Wronskian A p = — 2ip(l + p 2 ), independent of X. The zero of the Wronskian at 
p = i corresponds to the bound state mentioned above. Taking X > X', we obtain the 
Euclidean Green function discrete sum 



too 



%(X)¥(X' 



Gr " = 4^E- '- Wtf> (m). (18) 

Before proceeding, let us demonstrate that the Euclidean Green function is regular at 
the poles of the four sphere. This is a nontrivial check because the coordinates a and X 
are singular there, and the rescaling becomes divergent too. In the large X, X' limit, (|P8|) 
becomes 

00 -1 

Gf 3 \p, X, X') = 2k £ -e^x-x'lwffi'Qi) (19) 

n=3 U 



For n > 3 the Gaussian hypergeometric functions F(3 + n, 3 — n, 7/2, z) that constitute 

(n) 



the bitensor W,^\ l 3 have a series expansion that terminates, and they essentially reduce to 



Gegenbauer polynomials C^_ 3 (l — 2z). Using then the identity |T3[ 

00 _ 

^C^x)q l = (l-2xq + q 2 Y U (20) 
1=0 

with q = e _( ^ x ~ x '\ one easily sees that the sum ( |19"D indeed converges. 

We have the Euclidean Green function defined as an infinite sum (K). However, the 
eigenspace of the Laplacian on H 3 suggests that the Lorentzian Green function is most 
naturally expressed as an integral over real p. To do so we must extend the summand 
into the upper half p-plane. We have already defined the wavefunctions *$? P (X) as analytic 
functions for all complex p but we need to extend the bitensor as well. When the Green 
function is expressed as a discrete sum, it involves the bitensor W^ 1 3 (//) evaluated at p = ni 
with n integral. At these values of p, the bitensor is regular at both coincident and opposite 
points on S 3 , that is at p, = and p = tt. However, if we extend p into the complex plane 
we lose regularity at p — 0, essentially because the bitensor obeys the differential equation 
(O) with a delta function source at p, — 0. Similarly we must maintain regularity at p = n, 
since there is no delta function source there. The condition of regularity at 7r imposed by 
the differential equation for the Green function is sufficient to uniquely specify the analytic 
continuation of 3 (p) into the complex p-plane. The continuation is described in the 

Appendix, and the extended bitensor W^ 313 (p) is defined by equations (|A4|) and (|A7|) . 
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FIG. 1. Contour for the Euclidean Correlator. 



Now we are able to write the sum in (|T8| ) as an integral along a contour C\ encircling the 
points p = 3i, 4i, ...Ni, where N tends to infinity. For X > X' we have 



dp %(x)¥ p ( x') 

C\ psinhpn (1 + p 2 



A> A* 



G^ 3 (»,X,X f ) = K I *), \ „ p 2 \ ' W$> (n). (21) 



To see that (|2l|) is equivalent to the sum (|T8|) introduce 1 = coshj97r/ coshp7r into the 
integral. Then note that cothp7r has residue 7r _1 at every integer multiple of i. Finally, use 
( |A10|) to rewrite W^ tJ (fi) in the form regular at ji = used in (pT8|) . The factor of coshjwr 



from ( |A10| ) cancels that in the integrand. 



We now distort the contour for the p integral to run along the real p axis (Figure 1). 
At large imaginary p the integrand decays exponentially and the contribution vanishes in 
the limit of large N. However as we deform the contour towards the real axis we encounter 
two poles in the sinhp7r factor, the latter at p = i becoming a double pole due to the 
simple zero of the Wronskian. For the p = 2i pole, it follows from the normalisation of the 
tensor harmonics that W^ J = 0. Indirectly, this is a consequence of the fact that spin-2 
perturbations do not have a monopole or dipole component. At p = i we have a double 
pole, but although the relevant Schrodinger operator possesses a bound state, it does not 
generate a 'super-curvature mode'. Instead the relevant mode is a time-independent shift in 
the metric perturbation which may be gauged away We conclude that up to a term 

involving a pure gauge mode, we can deform the contour C\ into the contour shown in Figure 
1. For the moment, since the integrand involves a factor psinhp7r which has a double pole 
at p = 0, we leave the contour avoiding the origin on a small semicircle in the upper half 
]9-plane. 



Finally, in order to deal with the pole at p = 0, we re-express the integrand in fl2lD as 
a sum of its p-symmetric and p-antisymmetric parts. Denoting the integrand by I p we then 
have 

Gf 1 ' = \J dp(I p + J_ p ) + \ J dp(I p - J_ p ), (22) 

where the integral is taken from p = —oo to oo along a path avoiding the origin above. But 
/ dpl-p along this contour is equal to the integral of I p taken along a contour avoiding the 



origin below. The second term is therefore equal to the integral of I p along a contour around 
the origin. Hence we have 

- J dp(I p - J_ p ) = -7rzRes(J p ;p = 0). (23) 

We defer a detailed discussion of this term to Section V, because its interpretation is 
clearer in the Lorentzian region. Hence for the time being we just keep it, but it will turn 
out that it represents a non-physical contribution to the graviton propagator. 

In the p-symmetric part of the correlator, we can leave the integrand as a sum of I p and 
J_p. We henceforth denote the path from — oo to +00 avoiding the origin above by TZ. This 
shall turn out to be a regularised version of the integral over the real axis. Our final result 
for the Euclidean Green function then reads 

G%, M X, X>) = 1 1 V ^J^ + 

-iriRes{I p ;p = 0). (24) 



IV. TWO-POINT TENSOR CORRELATOR IN OPEN DE SITTER SPACE 

The analytic continuation into open de Sitter space is given by setting a = it and the 
polar angle Q = —ix- Without loss of generality we may take one of the two points, say Q' 
to be at the north pole of the three-sphere. Then fi = Q, and fi continues to — i\- We then 
obtain the correlator in open de Sitter space where one point has been chosen as the origin 
of the radial coordinate \. The conformal coordinate X continues to conformal time r as 
X = -r-f (see@). 

Hence the analytic continuation of the Euclidean mode functions is given by 

- -eftt^r) and *j,(X) - -e^^ p (r) (25) 
where the Lorentzian mode functions are 

tf£( r ) = (cothr + ip)e~ ipT . (26) 

They are solutions to the Lorentzian perturbation equation K^/p(r) = p 2 ^/^ (r). 

In order to perform the substitution fi = —ix, where x is the comoving separation on 
H 3 , we use the explicit formula given in the appendix for the bitensor regular at ji — 7T. 



The continued bitensor wMj,(x) is defined by the equations (|A7|), ( |A11| ) and (|A12| ). It 



can be seen from (|A12|) that it involves terms which behave as e ±p ^ x+7T \ One must extract 



the e p?r -factors in order for the bitensor to correspond to the usual sum of rank-two tensor 
harmonics on the real p-axis. To do so we use the following general identity. For t' — t > 0, 
we have (up to the p = i gauge mode) 

**» (T) *- (T 'W( P )=0, (27) 



c p (1 + p 1 



S 



where F(p) are the p-dependent coefficients occurring in the final (Lorentzian) form of the 
bitensor given in (| A 1 3|) . This identity follows from the analyticity of the integrand. By 
inserting 1 = sinhp7r/ sinhp7r under the integral, it is clear that the integral fl2"7|) with a 
factor e p7T / sinhp7r inserted equals that with a factor e~ pn / sinhp7r inserted. The resulting 
identity allows us to replace the factors e +p ^ x+7T ^ in the bitensor by e p ^ x ~^\ and vice versa 
in the analog integral of I- p closed in the lower half p-plane. 

For the tensor correlator we also need to restore the factor m _1 (r) to tij. It is convenient 
to define the eigenmodes $p(r) = SSfp (r)/a(r). The extra minus sign hereby introduced 
is cancelled by a change in sign of the normalisation factor Q p of the bitensor, which then 
becomes +(p 2 +4) / (307T 2 ). This corresponds to requiring the spacelike metric to have postive 
signature. We finally obtain the Lorentzian tensor Feynman (time-ordered) correlator, for 
t' - r > 0, 

-7rzRes(/^;p = 0), (28) 



where the Lorentzian bitensor W^M, is defined in the Appendix, equations (|A4j) and ( |A13|) 



In this section, we concentrate on the first term in (^), the integral over p, and ig- 
nore for the moment the second, discrete term. We first extract the symmetrised part, 
{{tij(x),ti'j'(x')}) , which is just the real part of the Feynman correlator. The imaginary 
part involves an integrand which is analytic for p — > 0: 



(U 3 (x),U>A*')) = \\ n ^ I ^2)^^'(x)cothp7r[^(r)$^(r') + $ L p (r)^(r')] 



-2k I dp 

Jo (1+p 2 



(29) 



It is straightforward to see that if we apply the Lorentzian version of the perturbation 
operator K to (|29| ) with an appropriate heaviside function of r — r', the imaginary term will 
produce the Wronskian of $^ p (r) and which is proportional to ip, times 5{r — r'). 

Then the integral over p produces a spatial delta function. From this one sees that our 
Feynman correlator obeys the correct second order partial differential equation, with a delta 
function source. The delta function source term in (^TJ) goes from being real in the Euclidean 
region to imaginary in the Lorentzian region because the factor ^Jg continues to iy/—g. 

The integral in (^) diverges as p~ 2 for p — > 0, in contrast with realistic models for 
inflationary universes where a reflection term in (|29| ) regularises the spectrum . However, 
as we immediately show, even in perfect de Sitter space the integral over p is perfectly 
finite. We rewrite the symmetrised correlator as an integral over real < p < oo as follows. 
Because the integrand in (|29|) is even in p, we have 



dp p7rcothp7r, 



2k 



^(r)^(r')W^> f ( X ) + 0(e), (30) 

ire 
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the second term being the contribution from the small semicircle around p = 0. Both terms 
may be combined under one integral. The resulting integrand is analytic as p — > and one 
can safely take the limit e — > 0. The symmetrised correlator is then given by 



(Vcoth^ 



7rp 2 V (1 + p 2 ) 



({tij(x),ti>y(x')}) = 
*f (T)^p(0] W^,(x) - 4i{(r)4}(0^(x)) , (31) 



where the Lorentzian bitensor W^fl is defined in the Appendix, equations (|A4j) and (|A13|) 



In this integral it may be written as 



W!$,(X) = E q^^^^r- (32) 

Vim 



The functions q£ (O) are the rank- two tensor eigenmodes with eigenvalues 
A p = — (p 2 + 3) of the Laplacian on if 3 . Here V = e,o labels the parity, and / and m are the 
usual quantum numbers on the two-sphere. At large p, the coefficient functions of the 
bitensor (see Appendix A) behave like psinpx- Hence the above integral converges at large 
p, for both timelike and spacelike separations. Furthermore, the correlations asymptotically 
decay for large separation of the two points. 

Equation (28), with the first term given by ( |5T] ) is our final result for the two-point 



tensor correlator in open de Sitter space, with Euclidean no boundary initial conditions. 
Contracting the propagator with the harmonics q^ elm and integrating over the three sphere 
reveals that the second term leaves the spectrum completely unchanged apart from cancelling 
the (divergent) contribution from the p 2 = divergence in the first term. We defer a detailed 
discussion of this result to the next section, in which we will also clarify the difficulties of 
the previous work on the graviton propagator in open de Sitter spacetime |fJ-Q . 

As an illustration let us compute the Sachs- Wolfe integral and show that all the mul- 



tipole moments are finite. The contribution of gravitational waves to the CMB anisotropy 
in perfect de Sitter space is given by 

0) = ~\ jH drt XXjT (r, X , 6, 0)| x=ro - T , (33) 

where r is the observing time. The temperature anisotropy on the sky is characterised by 
the two-point angular correlation function C(j), where 7 is the angle between two points 
located on the celestial sphere. It is customary to expand the correlation function in terms 
of Legendre polynomials as 

I ST ST \ 00 21 + 1 
C(l) = (^(0)^(7)^ = E ^C lPl (co S1 ). (34) 

Hence, inserting the Sachs- Wolfe integral into (0) and substituting (|3"ID for the two-point 
fluctuation correlator yields the multipole moments 



- u {wTp^-f www] 

-^o(t)^(t')QM'x') . (35) 
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In this expression we have written the normalised tensor harmonics qtyl elm (x, 0) as 
Qt x {x)Yirn{OA), where 



Ni(p) 



p 2 (p 2 + 1 



-(sinhx) 



1-2 



' -I d" 

sinhx d\ J 



i+i 



(cospx) 



(36) 



and 



Ni{p) 



(/-!)/(/ + !)(/ + 2) 
vrnU(j 2 + P 2 ) 



.1/2 



(37) 



It can readily be seen that the multipole moments are finite. With the aid of the explicit 
expressions and the wavefunctions (^) they can be numerically computed. 



V. CONCLUSIONS 

We have computed the spectrum of primordial gravitational waves predicted in open de 
Sitter space, according to Euclidean no boundary initial conditions. The Euclidean path 
integral unambiguously specifies the tensor fluctuations with no additional assumptions. 
The real space Euclidean correlator has been analytically continued into the Lorentzian 
region without Fourier decomposing it, and we obtained an infrared finite two-point tensor 
correlator in open de Sitter space, contrary to previous results in the literature @-@]. 

Let us now elaborate on the second, regularising term in the symmetrised correlator ( |3~T|) 
and the discrete p = contribution to the Feynman correlator given from the last term in 
(24). Not surprisingly, they have a similar interpretation. Their angular part W L i j%,(x) is 
equal to the sum of the tensor harmonics with eigenvalue X p (p = 0) = —3 of the Laplacian on 
H 3 . It has been known that a degeneracy appears between p 2 = tensor modes and pi = —4 
scalar harmonics [[J. More specifically, one has q^j = (VjVj — §7ijV 2 ) qi 2l ) lm where 
q(2i)im _ p^ lm Yi m . The discrete p 2 = tensor harmonics are the only transverse traceless 
tensor perturbations that can be constructed from a scalar quantity. But as a consequence 
of this, they are sensitive to scalar gauge transformations. Consider now the coordinate 
transformation £ a = (0, e$g (r)V l g < - 2 ^ /m ). Under this transformation the transverse traceless 
part of the metric perturbation in the perturbed line element (0) changes exactly by 
et\f lm = e$Q (r)q^ . Using the transverse-traceless properties of Uj it is easily seen that 
the action for tensor fluctuations is invariant under such transformations. Hence this tensor 
eigenmode is non-physical and can be gauged away. Note that since the functional form of £ is 
completely fixed this corresponds to a global transformation, analogous to the transformation 
— > 0+ constant for a massless field. To compute the Green function for a massless field 
one has to project out this homogeneous mode, and it is necessary to do the same here. One 
should therefore disregard the contribution from the discrete term in (24) to the Lorentzian 
correlator. This was actually also done in our computation of the tensor fluctuation spectrum 
about 0(4) instantons []IJ, although in that case not because the mode was pure gauge, but 
because it couples to the inflaton field, and is not represented by a simple action of the 
form (8). If a scalar field is present, the mode is most simply treated as a part of the scalar 
perturbations, as was done in 0. 
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In our result fl3~ID for the symmetrised correlator, the discrete gauge mode is set to zero 
because the second term cancels exactly the contribution from the p 2 = mode implicitly 
contained in the continuous spectrum. This automatic cancellation does not happen in the 
conventional mode-by-mode analysis where, if one chooses the most degenerate continuous 
representation of the isometry group 0(3, 1) of the hyperboloid if 3 , corresponding to the 
range p G [0, oo), one obtains a divergent correlator. 

It is clear that the underlying reason for these subtleties has to do with the different 
nature of tensor harmonics on compact and non-compact spaces. Hence, we could have 
expected the generation of the two discrete gauge modes simply from the analytic contin- 
uation of the completeness relation (14) of the harmonics on S 3 . Apart from the sum of 
the complete set of modes that constitute the delta function on H 3 , one obtains also three 
extra terms W%£\ 3 (aO, W)^ 3 (//) and W%\ 3 (//). The first term is zero, and the remaining 
two terms should respectively be viewed as sums of vector - and scalar harmonics. On the 
other hand, the fact that the scalar/tensor degeneracy appears precisely at the lower bound 
of the continuous spectrum is a peculiar feature of three dimensions. In the analogous com- 
putation in four dimensions for instance |IE |, this degeneracy happens at p 2 = —1/4 and 
consequently, there is no regularising term in the correlator. 

There is yet another way in which the exclusion of the degenerate modes from the 
perturbation spectrum can be interpreted. Remember that in non-compact spacetimes the 
decomposition (|6|) is only uniquely defined for bounded perturbations. Hence, the only 
way there can appear a degeneracy between the different types of fluctuations is for the 
degenerate modes to be unbounded. Indeed, on the three-hyperboloid the scalar p 2 = —4 
modes describe divergent fluctuations because the scalar spherical harmonics q( 2t ) lm grow 
exponentially with distance. The action of the above tensor operator renders only the 
^xi m com P onen ts of qlf lm finite at infinity. The remaining components still diverge as 



e x and correspond to exponentially growing fluctuations at large distances^. Since in 
cosmological perturbation theory one assumes the perturbation hij to be small, one must 
expand correlators in bounded harmonics. 

We want to emphasize that the regularity of the two-point tensor correlator does not de- 
pend on the Euclidean methods used in our work. One could have equally well computed the 
correlator on closed Cauchy surfaces for the de Sitter space where the subtleties encountered 
here do not arise, assuming the standard conformal vacuum for that slicing. One would then 
analytically continue the result to the open slicing. On the other hand, the Euclidean no 
boundary principle is an appealing prescription which avoids the arbitrary choice of vacuum 
otherwise needed. The path integral effectively defines its own initial conditions, yielding 
a unique and infrared finite Green function in the Lorentzian region. The initial quantum 
state of the perturbation modes, defined by the no boundary path integral, corresponds to 
the conformal vacuum in the Lorentzian spacetime. This is in many ways the most natural 
state in de Sitter space, but the regularity of the graviton propagator is independent of this 



The confusion arises because, due to the form of the metric inverse, scalar invariants are finite at 
infinity, e.g. qijq l] ~ e~ 2x . This also explains why the coefficient functions w^\x) m the bitensor 
Wij^jt asymptotically decay. 
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choice. The most important technical advantage of our method is that we deal throughout 
directly with the real space correlator, which makes the derivation independent of the gauge 
ambiguities involved in the mode decomposition. 

Finally, let us conclude by comparing the gravitational wave spectrum in perfect open 
de Sitter spacetime with the spectrum in realistic open inflationary universes. In both the 
Hawking-Turok and the Coleman-De Luccia model for open inflation there is an extra 
reflection term in the correlator because 0(5) symmetry is broken on the instanton [|TJ. 
This term gives rise to long-wavelength bubble wall fluctuations in the Lorentzian region. 
At first sight, the wall fluctuations seem to regularise the spectrum. However, adding and 



subtracting the second term in ([IT]) to the two-point tensor correlator in the 0(4) models (eq. 
(34) in and comparing that with our result ([H|) reveals that the wall fluctuations actually 
appear as an extra long-wavelength continuum contribution on top of the spectrum in perfect 
de Sitter space. Hence in both the Hawking-Turok and Coleman-De Luccia model there is 
an enhancement of the fluctuations compared to the perturbations in perfect de Sitter space. 
But the singularity in Hawking-Turok instantons suppresses the wall fluctuations because 
it enforces Dirichlet boundary conditions on the perturbation modes [I]]. Hence we expect 
the spectrum in perfect de Sitter space to be quite similar to the spectrum predicted by 
singular instantons. On the other hand, Coleman-De Luccia models typically predict large 
wall fluctuations, yielding a very different CMB anisotropy spectrum on large angular scales. 
The tensor fluctuation spectrum therefore potentially provides an observaional discriminant 



between different theories of open inflation |T5 |. 
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APPENDIX A: MAXIMALLY SYMMETRIC BITENSORS 



A maximally symmetric bitensor T is one for which a*T = for any isometry o of the 
maximally symmetric manifold. Any maximally symmetric bitensor may be expanded in 
terms of a complete set of 'fundamental' maximally symmetric bitensors with the correct 
index symmetries. For instance 



Tiji'j' 



h{lj)9ij9i'r + n i g ji ,n j , + n j g ii ,n j , + 71^--,% + n.g^n^ 



9ii'9jj> + 9ni'9; 



+ t 4 (/x)n i n i ri i ,n i , + t s (ji) g^n^, + n^^j, (Al) 



where the coefficient functions tj(n) depend only on the distance /x(f2, Q') along the shortest 
geodesic from Q to Q'. 71^(^,0') and 7^(0, ft') are unit tangent vectors to the geodesies 
joining Q and Q' and gij>(fl, fl') is the parallel propagator along the geodesic; V % g\ is the 



vector at £T obtained by parallel transport of V 1 along the geodesic from to Q,' fT7 



The set of tensor eigenmodes on S 3 or H 3 forms a representation of the symmetry group 
of the manifold. It follows in particular that their sum over the parity states V = {e, o} and 
the quantum numbers I and m on the two-sphere defines a maximally symmetric bitensor 
on S 3 (or H 3 ) [0] 



Viva 



(A2) 
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On S 3 the label p = 3i, Ai, ... It is related to the usual angular momentum k by p = i(k + 1). 
The ranges of the other labels is then < I < k and —l<m<l. On H 3 there is a 
continuum of eigenvalues p G [0, oo). We will assume from now that the eigenmodes on are 
normalised by the condition 



(A3) 



(p)Vlm 



The bitensor v ^ (//) appearing in our Green function has some additional properties 

arising from its construction in terms of the transverse and traceless tensor harmonics q. 

The tracelessness of wMa, allows one to eliminate two of the coefficient functions in ( |A1| ) 
It may then be written as 



(p) 



3m n a 



9i>j> ~ n v n v 



^n [i g jW n jl) +An i n j n i ,n j 



9a'9jj' + 9ji'9 



2n i g i , j ,n j - 2n if g ij n j , + Qn^^n^ 



(A4) 



This expression is traceless on either index pair ij or i'j'. The requirement that the bitensor 



be transverse V % Wjfhi 



and the eigenvalue condition (A 3 — \ P )W, 



(p) 



= impose addi- 
tional constraints on the remaining coefficient functions Wj P \fi). To solve these constraint 
equations it is convenient to introduce the new variables on S 3 (on H 3 , // is replaced by 
-ip) 



\ a(ji) 



7 da(fi) 
(p 2 +9) sin ft dfj, 



»/ 



In terms of a new argument z = cos 2 (/x/2) (or its continuation on H 3 ) 
and eigenvalue conditions imply for a(z) 



z(l 



. cPa(z) 



1-7, 
2 



da(z) 
dz 



{p 2 + 9)a(z) 



(A5) 

the transversality 
(A6) 



and then for the coefficient functions 



Wl = Q p ([2(\ p - Q)z{z - 1) - 2] a{z) + f [(A p + 6)z(z -\){z- 1)1 
w 2 = Q p [2{1 - z) {{X p -6)z + 3]a{z)-l[{X p + 6) z{z-l){z-'l))f3{z) 
w 3 = Q p ([-2(Ap - 6)z(z - 1) + 3] a{z) - | [(X p + 6)z(z - \){z - 1)] f3(z)) 



(A7) 



with A p = (p 2 + 3). 



The above conditions leave the overall normalisation of the bitensor undetermined. To 
fix the normalisation constant Q p we contract the indices in the coincident limit z — > 1. 
This yields |L8| 



(p) ij 



•J 



J2 q^ )Vlm iSl)q^ Vlm 



30Q p a(l). 



Vim 



By integrating over the three-sphere and using the normalisation condition 
tensor harmonics one obtains Q p = — 30 ^2^i) • 



(A8) 
on the 
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Notice that (|A£]) is precisely the hypergeometric differential equation, which has a pair 
of independent solutions a(z) = 2^1 (3 + ip, 3 — ip, 7/2, z) and 

a(l — z) = 2-^i(3 + ip,3 — ip, 7/2, 1 — z). The former of these solutions is singular at z = 1, 
i.e. for coincident points on the three-sphere, and the latter is singular for opposite points. 
The solution for f3(z) follows from (|A5|) and is given by 



P(z)= 2^(4-^,4 + ^,9/2,2). (A9) 
The hypergeometric functions are related by the transformation formula (eq. [15.3.6] in |fL9| ) 

r(c)r(c — a — b) . . 

2 Fi(a, 6, c, z) = — — 2 Fi(a, b, a + b - c, 1 - z) 

T(c — a)T(c — b) 

+ T{C) rf\rL~ C) ^ ~ z y~ a ~ b 2F 1 (c-a,c-b,c-a-b,l-z). (A10) 



Only for the eigenvalues of the Laplacian on S 3 , i.e. p = in (n > 3), the term on the 
second line vanishes for 2-^i(3 + ip, 3 — ip, 7/2, z). For these special values, a(z) and a{l — z) 
are no longer linearly independent but related by a factor of (— l) n+1 , and they are both 
regular for any angle on the three- sphere. In fact, the hypergeometric series terminates for 
these parameter values and the hypergeometric functions reduce to Gegenbauer polynomials 
Cnl 3 (l — 2z). We have a choice between using a(z) and a(l — z) in the bitensor for these 
values of p. Since F(l — z) — > 1 for coincident points, it is more natural to choose a(l — z) 
in the bitensor appearing in the Euclidean Green function However, to obtain the 

Lorentzian correlator, we had to express the discrete sum ([18]) as a contour integral. Since 
the Euclidean correlator obeys a differential equation with a delta function source at fi = 0, 
we must maintain regularity of the integrand at fi = 7r when extending the bitensor in the 
complex p-plane. In other words, for generic p, we need to work with the solution a{z), 
rather than a(l — z). Therefore, in order to write the Euclidean correlator as a contour 
integral, we first have replaced F(l — z) by F(z)(— l) n+1 , by applying ( A10|) to ( fL8|) , and we 



then have continued the latter term to — (coshjwr) 2 F\ (3 + ip, 3 — ip, \,z). 

We conclude that the properties of the bitensor appearing in the tensor correlator com- 
pletely determine its form. Notice that in terms of the label p we have obtained a 'unified' 
functional description of the bitensor 3 on S 3 and H 3 . Its explicit form is very different 
in both cases however, because the label p takes on different values. But it is precisely 
this description that has enabled us in Section IV to analytically continue the correlator 
from the Euclidean instanton into open de Sitter space without Fourier decomposing it. We 
shall conclude this Appendix by giving the explicit formulae for the coefficient functions of 
the bitensor W*£fj, appearing in our final result (|3l|) . With this description, they can be 
obtained by analytic continuation from S* 3 . 

To perform the continuation to if 3 we note that the geodesic separation [i on S 3 continues 
to —i\ where \ is the comoving separation on H 3 . Hence the hypergeometric functions on 
H 3 are defined by analytic continuation (eq. 15.3.7 in [19|]) and may be expressed in terms 
of associated Legendre functions as 

(AH) 




(- sinh x)- 5/2 P:$ +ip (- cosh X ), 
15 J% (- sinh x)- ll2 PZl'/ 2 2 +ip {- cosh x). 
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Using the relation — cosh(x) = cosh(x — m), the Legendre functions on H 3 may be expressed 
as 



-5/2 , 
-l/2+ip\ 



-7/2 , 
-l/2+ip\ 



cosh x) 



cosh x) 



^T— 1 +P 

— 7rsmnx v 

i siul. ,•[ i y;+7r) 

p 2 



2p 



-\ 1 (4 + j> 2 ) 1 [—3 cothx cosh p{n + ix) 
p 2 )(l + coth 2 x) + (4 + p 2 )cosech 2 x 



-(1 



coshp(7r + ix)(p 2 — 11 — 15cosech 2 x) 



x 



(A12) 



-6 



i sinhp(ix+ 



p 2 ) coth 3 x + (p 2 + |) cothx cosech 2 ^ 



In the text, we have extracted the factors e ±p7T in these expressions in order to make contact 
with the usual description of the tensor correlator in terms of tensor harmonics on if 3 . The 
coefficient functions of the bitensor W L M jt (%) in our final result (j3~ID for the tensor correlator 



arc 



w 2 



w 3 



cosech 5 x 

47T 2 (p 2 +l) 



V 



(3 + (p 2 + 4) sinh 2 x - P 2 {p 2 + 1) sinh 4 x) 



cospx(3/2 + (p 2 + 1) sinh x) sinh 2% 



cosech 5 x 

4vr 2 (p 2 +l) 



smpx 
V 



(3 + 12 cosh x ~ 3p 2 (l + 2 cosh x) sinh 2 x 



cosech ° x 
4vr 2 (p 2 + l) 



-p > (p + 1) sinh x) + cos PX(— 12 — 3 cosh x 
+2(p 2 — 2) sinh 2 x + 2(p 2 + 1) cosh x sinh 2 %) sinh x 

™ (3 - 3p 2 sinh 2 x + P 2 {P 2 + 1) sinh 4 x) 
+ cospx( — 3/2 + (p 2 + 1) sinh 2 x) sinh 2% 



(A13) 



As mentioned before, the bitensor W^-fj, equals the sum (|A2|) of the rank-two tensor eigen- 
modes with eigenvalue \ = —(p 2 + 3) of the Laplacian on H 3 . For x ~ * these functions 
converge and they exponentially decay at large geodesic distances. 
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